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ABSTRACT 

We have conducted a search for giant pulses from four millisecond pulsars 
using the 100 m Green Bank Telescope. Coherently dedispersed time-series from 
PSR J0218-I-4232 were found to contain giant pulses of very short intrinsic du¬ 
ration whose energies follow power-law statistics. The giant pulses are in phase 
with the two minima of the radio integrated pulse prohle but are phase aligned 
with the peaks of the X-ray prohle. Historically, individual pulses more than 10- 
20 times the mean pulse energy have been deemed to be “giant pulses”. As only 
4 of the 155 pulses had energies greater than 10 times the mean pulse-energy, 
we argue the emission mechanism responsible for giant pulses should instead 
be dehned through: (a) intrinsic timescales of microsecond or nanosecond dura¬ 
tion; (b) power-law energy statistics; and (c) emission occurring in narrow phase- 
windows coincident with the phase windows of non-thermal X-ray emission. Four 
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short-duration pulses with giant-pulse characteristics were also observed from 
PSR B1957-I-20. As the inferred magnetic helds at the light cylinders of the 
millisecond pulsars that emit giant pulses are all very high, this parameter has 
previously been considered to be an indicator of giant pulse emissivity. However, 
the frequency of giant pulse emission from PSR B1957-I-20 is signihcantly lower 
than for other millisecond pulsars that have similar magnetic fields at their light 
cylinders. This suggests that the inferred magnetic held at the light cylinder is 
a poor indicator of the rate of emission of giant pulses. 

Subject headings: pulsars:general — pulsarsdndividual (PSR J0218-I-4232, PSR J1012-I-5307, 
PSR J1843-1113, PSR B1957+20) 


1. Introduction 

The Crab radio pulsar was discovered through the direct detection of strong individual 
pulses (Staelin & Reifenstein 1968). Further studies revealed that the strongest pulses fol¬ 
lowed power-law energy statistics (Argyle & Gower 1972) distinct from the Gaussian statistics 
of the general pulse population (Cordes 1976). In an observation by Lundgren et al. (1995) 
around one in 1200 pulses had an energy greater than 20 times the mean pulse energy, (E). 
Despite this, Cordes et al. (2004) found that at all radio frequencies phase-coherent summa¬ 
tion of the giant pulses gives a higher signal-to-noise ratio than summation of all the pulses. 
Extraordinarily, the giant pulses also have structure that is signihcantly narrower than the 
mean pulse. Hankins et al. (2003) observed pulses with that had structure persisting for less 
than 2 ns and inferred that the brightness temperatures of these pulses are Tb lO^^K. 

The young Crab-like pulsar B0540—69 in the Large Magellanic Cloud also emits gi¬ 
ant pulses (Johnston & Romani 2003). In 31.2hr of observations at a center frequency of 
1390 MHz Johnston et al. (2004) only detected the integrated emission-prohle of PSR B0540—69 
at a very low level of signihcance. Despite their difficulty in detecting the integrated emis¬ 
sion, Johnston et al. were able to detect and analyse 141 individual pulses. The relative ease 
with which giant pulses can be seen over large distances has led several authors to advocate 
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their detection as a way to find extra-galactic pulsars (see e.g. Johnston & Romani 2003; 
Cordes et ah 2004). 

To date, no other young pulsars have been found to emit pulses with the high energies 
and extremely short durations characteristic of the giant pulses from the Crab pulsar. Three 
young pulsars have been found to emit narrow pulses of emission showing power-law statistics 
(Johnston et al. 2001; Johnston & Romani 2002; Cairns et ah 2004). However, it is not clear 
that the pulses should be classed as true “giant pulses” because the power-law tails have 
only been seen to extend to low energies. In addition, the structure of these events has thus 
far not been shown to have timescales as short as those of giant pulses from the Crab pulsar. 

The recycled pulsars B1937-I-21, B1821—24, and J1823—3021A also emit giant pulses 
despite having markedly different periods (P) and period derivatives (P ) to the Crab pulsar 
(Cognard et al. 1996; Romani & Johnston 2001; Knight et al. 2005). One common factor 
between these millisecond pulsars , PSR B0540—69, and the Crab pulsar is that they all have 
very high magnetic fields inferred at their light cylinders^ Plc oc p-2-5po.5 viewed 

in the context of the known millisecond pulsar population, the three giant pulse emitters 
also have very low characteristic ages r = P/ (2P) and very high spin-down luminosities 
E oc P~^P. PSRs B1821—24, B1937-I-21, and J0218-I-4232 have some of the highest X- 
ray luminosities of all millisecond pulsars (Becker & Triimper 1999; Grindlay et al. 2002; 
Cusumano 2004; Heinke et al. 2005). The emission from all three pulsars is non-thermal, 
and the X-ray profiles of PSRs B1821—24 and B1937-I-21 align in phase with their giant 
pulse emission (Romani & Johnston 2001; Cusumano et al. 2003). Another field pulsar with 
a high X-ray luminosity is PSR B1957-I-20 (Becker & Triimper 1999). However, no X-ray 
pulsations have been detected from this source and it is unclear how much of the emission 
originates from the bow-shock between the pulsar wind and the companion wind (Stappers 
et al. 2003). 

In this paper we present the results of a sensitive baseband search for microsecond- 
timescale emission from four millisecond pulsars. Upper limits are placed on emission 
from PSRs J1843—1113 and J1012-I-5307, and a new population of short-duration pulses 
is reported for PSR B1957-I-20. A previously unknown population of giant pulses from 
PSR J0218-I-4232 is characterized and the results used to clarify the defining characteristics 
of giant pulse phenomenology. 


■^The ATNF Pulsar Catalogue has been used to obtain the pulsar parameters and statistics used in this 
paper. See: http://www.atnf.csiro.au/research/pulsar/psrcat. 



2. Observations and Data Analysis 


All observations were taken using the 100 m NRAO Green Bank Telescope (GBT) from 
2004 August to 2005 January at frequencies in the ranges of 793-921 and 1341-1469 MHz. 
Data were acquired using the Galtech-Green Bank-Swinburne Recorder II (see Jacoby 2005). 
This instrument real-samples one or two dual polarization 64 MHz-wide bands at the Nyquist 
rate. Software algorithms similar to those described by van Straten (2003) were used to 
synthesize hlter banks. The first step of the technique is to Fourier transform the raw 
voltages to the frequency domain and divide the spectra into a series of sub-bands. Each 
sub-band is multiplied by an inverse-response hlter (kernel) for the interstellar medium (ISM) 
(see, e.g. Hankins & Rickett 1975; Stairs 1998). The sub-bands are then individually Fourier 
transformed back to the time domain to give numerous time series, each having coarser time 
resolution than the original. This technique avoids the spectral leakage suffered by forming 
the hlter bank hrst and then transforming back to the Fourier domain to deconvolve. By 
splitting the input signal into sub-bands the dispersive smearing that has to be accounted for 
is essentially reduced to that of an individual sub-band. This means that to hrst order the 
number of samples required for the initial forward transform is inversely proportional to the 
number of sub-bands in the hlterbank. Gonsequently forming such a “coherent hlterbank” 
uses much shorter transforms than single-channel coherent dedispersion. In practical terms 
this means the algorithm can use high-speed memory more exclusively and therefore is 
computationally faster. 

Goherent dedispersion and channel summing were repeatedly applied to cover a range 
of dispersion measures (DMs) typically within ±0.1pccm“^ of the published pulsar DM. 
This guaranteed that our sensitivity would never be reduced due to DM error. Data were 
square-law detected and combined to give a dataset with bandwidth of 64 or 128 MHz. 
These time series were then searched for broad-band emission by summing adjacent samples 
at time resolutions between 1 and 128 /rs. Any two samples with total flux ISa (Her for 
PSR J0218-I-4232) or more above the local mean were further reduced to produce candidate 
plots for human scrutiny. 

Table 1 summarizes the observations taken. Golumns 1-3 show the pulsar name, center 
frequency, and bandwidth respectively. Golumns 4 shows the observation duration, and 
column 5 shows the number of pulses observed. The mean pulse energy and 1/is sensitivity 
threshold are shown in columns 6 and 7 respectively. The last column shows the number of 
individual pulses detected. For PSR J0218-I-4232 this column shows the number of pulses 
detected in each of the “A” and “B” pulse-phase regions discussed in Section 3.1.1 and shown 
in Figure 1. The system equivalent flux densities for the frequency bands centered at 825-889 
and 1373-1437 MHz ranged between 12-14 and 9.1-9.4 Jy respectively. 
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3. Search Results 
3.1. PSR J0218+4232 

3.1.1. Properties of the Pulses 

A total of 155 emission events were detected from PSR J0218+4232. As these aligned 
in two distinct pulse-phase windows (see Figure 1) they are all identihed as individual pulses 
from PSR J0218-I-4232. Figure 1 also shows the phases of the pulses relative to X-ray 
prohle and the full-width, half-maximum (FWHM) energies of the pulses relative to the 
mean pulsed flux density. PSR J0218-I-4232 has a signihcant ~ 50% unpulsed component 
(Navarro et ah 1995) which was not accounted for in the calculation of the average pulse 
energy. X-ray emission is more prevalent in the earlier “A” phase window, but more giant 
pulses were detected at the later “B” phase window. Our observations therefore show that 
although giant pulses in the radio band appear to originate in the same part of the pulsar 
magnetosphere as X-ray emission, they are modulated by different processes. For the August 
857MHz observation the “A” and “B” emission windows spanned 81/xs (0.035 periods) and 
123 fj,s (0.053 periods) respectively. Similar widths of 3% and 4% of phase were measured for 
the windows at 1373 MHz. The phase regions in which PSR B1937-I-21 emits giant pulses 
are much narrower. At 1650 MHz its two windows are 10.7 and 8.2/is wide, or 0.007 and 
0.005 periods wide (Soglasnov et ah 2004). The giant pulses found on the main emission 
component of PSR J1823—3021A at 685 MHz have a similar phase range to the pulses of 
PSR J0218-I-4232 of about 0.04 periods or 220/xs (Knight et al. 2005). 

The Crab pulsar and PSR B1937-I-21 emit 10-20(77) pulses at high rates, and so energy 
thresholds in this range have been used to distinguish giant pulses from ordinary emission 
(Argyle & Gower 1972; Cognard et al. 1996). Only 3 of the 139 pulses seen at 857MHz 
from PSR J0218-I-4232 had energies greater than 10(77), and so this pulse population is 
not particularly strong compared to the giant pulse populations of the Crab pulsar and 
PSR B1937-I-21. However, the argument that these pulses arise from the same giant pulse 
emission mechanism is compelling. Firstly, the cumulative distribution of pulse energies 
shown in Figure 2 shows that the strongest pulses have power-law statistics. The tapering 
off at low energies is due to the widths of the pulses being underestimated due to noise. 
The pulses are very narrow and align in phase with the non-thermal X-ray pulses. All these 
properties are shared by the giant pulses of the Crab pulsar and PSR B1937-I-21. In addition, 
the pulses from PSR J0218-I-4232 occur at the minima of the integrated emission prohle and 
therefore do not contribute to the main emission components. Consequently they cannot 
be interpreted as strong “ordinary” pulses. The pulses from PSR J0218-I-4232 therefore 
demonstrate that the giant pulse phenomenon can no longer be dehned through arbitrary 
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bounds on pulse energy. Better phenomenological criteria are narrow pulse-widths, power- 
law statistics, and emission occuring in narrow phase-windows that align with non-thermal 
X-ray emission. 

Johnston & Romani (2004) also suggested that power-law statistics and emission at 
special phases were the dehning characteristics of giant pulses. Their hlter bank observations 
in previous work (Romani & Johnston 2001) were unable to constrain the width of the giant 
pulses from PSR B1821—24. Our observations have shown that PSRs J1823—3021A (see 
Knight et ah 2005) and J0218-I-4232 have intrinsically narrow pulses. We argue that giant 
pulses always have narrow widths, and that this property can be added to those presented 
by Johnston & Romani in dehning the giant pulse phenomenon. 

The fraction of pulses detected at phase “A” almost halved from 0.24 in August to 
0.13 in October. As the search did not discriminate on the basis of phase, this difference 
could be interpreted as being due to variation in the rate of giant pulse emission for each 
phase window. However, the rate change only becomes readily apparent when viewed in 
terms of the detection counts regardless of pulsar hux (see left panel of 2), and not when 
viewed in terms of energy relative to the mean pulse-energy (see right panel of 2). Small 
number statistics are therefore a more likely cause of the disparity- the non-detection of ~ 6 
low-energy pulses can explain the rate change. 


3.1.2. Comparison of Emission Rates 

The probability of a pulse having energy greater than Eq can be expressed as: 

P{E > Eo) = KEf^. (1) 

Here Eq is in units of the mean pulse energy. Integrating gives an expression for the fraction 
of pulse flux emitted in the form of giant pulses of energies greater than Eq\ 

SgAE > Eo) = -^El-P (2) 

The best hts for the 857 MHz pulses with energies greater than 25JyjU.s are shown in 
Table 2. No satisfactory £t was obtained for the October “A” pulses. Estimates of the 
relative rate at 1373 MHz and rates for other pulsars are also shown. The first three columns 
show the pulsar, center frequency, and the phase range the power law is valid for. Columns 4 
and 5 show the best fits for K and a. The probability that a pulse has E > 2D{E) is shown 
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in column 6. Columns 7 and 8 show the fraction of flux that is emitted as giant pulses of 
energies greater than 20(i7) and 0.1(77) respectively. 

The power-law energy distributions of the Crab pulsar, PSR B0540—69, and PSR B1937-I-21 
(at 430MHz) do not extend to energies as low as 0.1(77). Soglasnov et al. (2004) hnd that 
at 1650MHz the giant pulses from PSR B1937-I-21 extend to energies of 0.016-0.032 (77). 
The power-law exponents for the millisecond pulsars PSR B1821—24 and PSR J1823—3021A 
are poorly known. However, at 1400-1500 MHz they emit a giant pulse of more than 28(77) 
at frequencies of ~ 8.5 x 10“^ and ~ 4.6 x 10“® respectively (Romani & Johnston 2001; 
Knight et al. 2005). For comparison, the work of Soglasnov et al. gives for PSR B1937-I-21 
an emission rate of P{E > 28(77)) = 2.6 x 10“®. The observed pulse energy distribution 
is the product of the intrinsic distribution and the spectra of propagation effects such as 
interstellar scintillation. Scintillation is particularly strong for PSR B1937-I-21 on timescales 
of minutes at frequencies in the vicinity of 1-2 GHz and could potentially lead to different 
studies obtaining different results. Kinkhabwala & Thorsett (2000) obtain parameters for 
PSR B1937-I-21 at 1420MHz of a = 1.8 and P(77 > 28(77)) ~ 4.0 x 10“^ which are quite 
different from those found by Soglasnov et al. 

It is apparent in Table 2 that PSR J0218-I-4232 has a much lower rate of giant pulse 
emission than other giant pulse emitters. The total fraction of its pulsed energy emitted in 
the form of giant pulses with energies greater than 0.1(77) is about 0.1%. Such giant pulses 
can occur at rates of up to one per ~ 200 pulsar rotations. If the cut-off point of the power 
law occurs at 0.1(77) then the ~ 10% of the pulse prohle where the giants occur should have 
a flux enhancement caused by the giant pulses of ~ 1% of the mean flux density. If the power 
law extends to 0.01(77) then the flux enhancement increases to the 10% level. Extension to 
energies much lower than 0.01(77) does not seem plausible given the lack of large components 
in the emission regions of the giant pulses. 

The power-law £t for the ten most energetic pulses seen at 1373 MHz is summarized in 
Table 2. A 20(77) pulse at 1373 MHz is emitted about 1.7 times more frequently than the 
August 857 MHz “B” pulses. It should be noted that the formal uncertainty on a of ±0.1 
makes this estimate somewhat uncertain. 


3.1.3. Pulse Durations 

All the 857 MHz pulses had FWHM durations of < 3.2/rs. The strongest pulse had 
a FWHM duration of 2.6 fj,s. To investigate the possibility of substructure, this pulse was 
coherently dedispersed at a time resolution of 15.625 ns. The initial portion of the pulse is 



shown in the top panel of Figure 3. The hnite rise time is only resolved at sampling intervals 
less than 125 ns and persists if the DM is slightly altered. At high time resolution the 
noise statistics are better modelled using chi-squared distributions than with the standard 
Gaussian approximation. The noise statistics therefore become more positively skewed at 
higher time resolutions and so much of the substructure seen at 15.625 ns time resolution is 
likely to be spurious. 

The strongest pulse seen at 1373 MHz as shown in the middle panel of Figure 3 is 
signihcantly narrower. At a time resolution of 125 ns it is of order 500 ns wide. 

Strong spikes following the main emission peak persist for about 8 times longer at 
857 MHz than at 1373 MHz, and so the pulse widths are roughly consistent with the 
scaling law of Kolmogorov-spectrum interstellar scattering (Bhat et al. 2003). We think the 
hnite rise-time seen at both 857 and 1373 MHz is a consequence of propagation through a 
thick scattering screen (see, e.g. Williamson 1973) rather than intrinsic substructure. 


3.1.4- Timing of Giant Pulses 

The giant pulse emission of PSR J0218-I-4232 occurs over much narrower ranges of pulse 
phase than the integrated pulses. It is therefore important to consider whether timing of 
PSR J0218-I-4232 can be improved by timing only the giant pulses. We formed a standard 
prohle from the brightest giant pulse and cross correlated the giant pulses with it to obtain 
an arrival time. The 56 giant pulses in the August observation in phase range “B” that had 
arrival time errors less than 0.5 ps had an rms residual of 24ps. The error in arrival time 
for the whole group was therefore about 3 ps. However our timing of the mean prohle for 
this observation obtains an rms residual of 6 ps using 16.8 s integrations, which we would 
expect to improve signihcantly with increased integration. Therefore conventional timing 
gives superior results to timing using giant pulses. 


3.2. PSR J1012+5307 

PSR J1012-I-5307 is a 5.3 ms pulsar with a characteristic age of 8.6 Gyr (Lange et al. 
2001). It has a Rlc 68 times smaller than that of PSR B1937-I-21. No individual pulses 
were detected from PSR J1012-I-5307. Our result is consistent with the fact that to date no 
millisecond pulsars with low values of Rlc and large characteristic ages have been observed 
to emit giant pulses. Edwards & Stappers (2003) observed PSR J1012-I-5307 for 1800 s at 
1380 MHz using the Westerbork Synthesis Radio Telescope. With a sampling interval of 
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51.2/is they detected 70 individual pulses with energies of up to hve times the mean pulse 
energy. Our observations establish that it is very unlikely that the pulses uncovered by 
Edwards & Stappers have the short < 1 /xs timescales characteristic of the giant pulses of 
PSR B1937+21. 


3.3. PSR B1957+20 

PSR B1957+20 has the third highest Rlc of all millisecond pulsars and was therefore 
targeted by Knight et al. (2005) as a potential source of giant pulse emission. Knight et al. 
failed to detect any pulses in 7700 s of observations at a center frequency of 685 MHz using 
the Parkes Radio Telescope. It is well known that the pulsar wind of PSR B1957+20 causes 
gas to be ablated from its companion (Fruchter et al. 1990; Krolik & Sincell 1990). This 
ionized gas causes eclipses at orbital phases (0) near 0.25. Knight et al. had suggested that 
the gas could scatter-broaden any giant pulses beyond reasonable detection levels. However, 
signihcant broadening cannot occur at all orbital phases, as in 8003 s of observations using 
the GBT we detected four narrow pulses from PSR B1957+20. Our observations spanned 
0.40 < 0 < 0.67; the earliest pulse arrived at 0 = 0.41. 

To estimate the energies of the pulses we formed a 512-bin prohle of each pulse and 
calculated the FWHM energy. The pulses had energies of 4.5-8.6(77). At this coarse time 
resolution virtually all of the pulse flux for these pulses is encompassed in our estimate. 
Adjustment of the DM used for coherent dedispersion and channel summing causes changes 
in the noise characteristics of the on-pulse region. Peak intensity, pulse morphology, and 
pulse width all vary with DM, and so determination of the true DM and therefore true pulse 
shape becomes dependent on the exact criteria used to optimize DM. The bottom panel 
of Figure 3 shows the strongest pulse. Although the main pulse component appears very 
narrow, there is a very weak underlying emission region about it of microsecond duration. 
Other pulses optimize at DMs that differ by 0(10“^) pccm“^. We think that this DM 
uncertainty is caused by the low signal strength of the pulses. It means that the prohles 
shown at high time resolution do not necessarily represent the true pulse form. For weak 
pulses like these, we suggest that the true nature of the pulses in terms of their individual 
DMs, intrinsic widths, and substructure requires the DM to be accurately determined via 
multi-frequency observations. 

The four pulses fall in a narrow 20 /xs pulse window that covers the peak of the main emis¬ 
sion component (see Figure 4). The other pulse components of the integrated prohle are 50% 
or more weaker than the main component. It is reasonable to suppose that PSR B1957-I-20 
might emit narrow pulses similar to those observed that are phase-aligned with the other 
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components. If these pulses exist and are amplitude modulated in a similar fashion to the 
“ordinary” pulse emission, they would be 50% or more weaker than the main-component 
pulses we see. As our initial detection threshold was Ida and the main-component pulses 
were detected at 14-16(T, our observations do not place good bounds on the existence of such 
pulses. 


3.4. PSR J1843-1113 

PSR J1843—1113 is a solitary 1.8 ms pulsar with a characteristic age of ~ 3Gyr. Its Rlc 
is very high — about 0.2 times that of PSR B1937-I-21. No spikes of broad-band emission 
were detected from it, suggesting that if it does emit giant pulses they are very weak and/or 
infrequent. Because PSR J1843—1113 is close to the plane of the Galaxy and has a relatively 
high DM we cannot rule out that it emits pulses similar to those of PSR B1957-I-20, but 
which are scatter-broadened beyond our sensitivity limits. 


4. Discussion 
4.1. Pulse Populations 

Joshi et ah (2004) reported the detection of an unresolved ~ 129(77) (925 Jy/xs) large- 
amplitude pulse from PSR B1957-I-20 in observations centered at 610 MHz. As this pulse 
is much larger than any detected in our observations it is instructive to consider whether it 
constitutes: (a) a member of a different pulse population of longer intrinsic duration; (b) the 
very high-energy tail of the distribution we observed; or (c) some sort of noise event that is 
unrelated to the pulsar. If the pulse reported by Joshi et al. is as broad as their sampling 
time (258 /xs) then our 825 MHz detection threshold for summing two 128 /xs samples of 45(77) 
should have found similar pulses, and so our observations are not consistent with hypothesis 
(a). Alternatively if we assume the pulses follow a a = —1.4 energy distribution then our 
observations imply that a 129(77) pulse should be emitted on average once every 61 hr. 
Given that Joshi et al. observed PSR B1957-I-20 for < 1 hr and much steeper power-law 
distributions have been found for other giant pulse emitters (see, e.g. Romani & Johnston 
2001) we hnd hypothesis (b) untenable. The detection criterion used by Joshi et al. was 
that a pulse must exceed 3.5a in two bands. With an rms of ~ 1 Jy over their 16 MHz band 
it is apparent that a bare detection corresponds to 1300Jy/is, or 180(77). For a normally 
distributed noise floor approximately eleven noise spikes would be expected to exceed this 
threshold in their sample of ~ 10® pulses. The fact that multiple noise spikes with higher 
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energies than the pulse are expected to be present in the Joshi et ah data makes it difficult 
to argue that the pulse is not background noise. This in turn implies that we have presented 
the hrst evidence for a population of giant pulses from PSR B1957+20. 

Joshi et ah also reported the detection at 610 MHz of three unresolved large-amplitude 
~ 258 ps wide pulses from PSR J0218-I-4232 with energies of 48-51 {E). The event rate of 
this pulse population is P{E > 48(i7)) = 1.4 x 10“®, which is 40 times higher than our 
August rate for the “B” phase range of P{E > 48(77)) = 3.6 x 10“®. Our August detection 
threshold for summing two 128 ps samples of 7(77) means that we should have easily detected 
the Joshi et ah pulses. Therefore the pulses of Joshi et ah are not a separate pulse population 
that is simply stronger than the one we observed. Furthermore, the Joshi et al. pulses occur 
at a different phases to the pulses we saw, so the hypothesis that Joshi et al. were extremely 
fortunate in detecting the high-energy tail of our population is not at all plausible. Pulses 
similar to those reported by Joshi et al. should also have been seen by Edwards & Stappers 
(2003), who did not detect any pulses above 26(77) in an 1800s observation centered at 
328 MHz. If the noise floor of Joshi et al. is normally distributed then approximately 36 
noise spikes in their sample would be expected to exceed their criterion of 3.5(T in both 
bands and therefore have a similar energy to the pulses reported. The three pulses are 
therefore not distinguishable from background noise and are likely to be spurious. The 
only type of strong pulses not ruled out by our data reduction are those with timescales 
comparable to PSR J0218-|-4232’s 2.3 ms pulse period. However, such pulses probably would 
have had substructure detectable in our searches. It is more likely PSR J0218-I-4232 only 
emits one population of strong pulses and that is the population of giant pulses unveiled by 
our observations. 


4.2. Giant Pulses from PSR B1957-|-20 

The pulses seen from PSR B1957-I-20 have sub-microsecond timescales and are several 
times stronger than the mean pulse. All four coincide with the main emission component in 
a similar fashion to the giant pulses of PSR J1823—3021A. Even without evidence for power- 
law statistics it is tempting to categorize PSR B1957-I-20 as a giant pulse emitter. However, 
the pulses can also be explained as strong pulses of ordinary emission that are exceptionally 
narrow. This hypothesis is supported by the fact that pulses from PSR J0437—4715 exhibit 
an anti-correlation between pulse width and pulse strength (Jenet et al. 1998). The strongest 
“ordinary” pulses from PSR J0437—4715 are then much more readily detected in single¬ 
pulse searches, and therefore could potentially masquerade as a giant-like population. Jenet 
et al. found pulses as short as 10 ps in their ~ 3000 s of observations. Therefore it is not 
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unreasonable to suggest that in ~ 8000 s PSR B1957+20 could emit several “ordinary” pulses 
consisting of very short spikes superimposed on microsecond-timescale emission bursts. The 
fact that we did not detect any microsecond-timescale emission from PSR J1012-I-5307 means 
the pulse substructure seen by Edwards & Stappers (2003) is broader than that seen for 
PSR B1957-I-20. Similarly, Knight et ah (2005) did not hnd any substructure in pulses from 
PSR J1603—7202 as short as their 4 /j.s sampling time. Microstructure within ordinary pulses 
from millisecond pulsars therefore does not seem to have characteristic timescales as short 
as those of the PSR B1957-I-20 pulses. 

Insight into whether or not the pulses from PSR B1957-I-20 are plausibly “giant” can be 
gained by comparing the properties of PSR B1957-I-20 and pulsars that emit giant pulses. 
Table 3 summarizes the attributes of the millisecond pulsars previously known to emit giant 
pulses (top) and the pulsars we observed (bottom). Each of these two groups is sorted by 
right ascension. The first three columns show the pulsar name, period, and period deriva¬ 
tive respectively. The period derivatives have been corrected for kinematic effects where 
possible (Shklovskii 1970; Damour & Taylor 1991). PSRs B1821—24 and J1823—3021A are 
located within globular clusters, and so acceleration in the cluster potential will contribute 
to the observed P for these pulsars. The magnitude of the cluster contribution to P is very 
uncertain, but has been estimated to be < 0.06P for PSR B1821—24 (Phinney 1993) and 
< 0.7P for PSR J1823—3021A (Stappers 1997). The proper motions of PSRs J0218-I-4232 
and J1843—1113 are unknown, but a 100kms“^ velocity equates to a Shklovskii-term contri¬ 
bution to P of just 0.6% for PSR J0218-I-4232 and 10% for PSR J1843—1113. Columns 4-7 of 
Table 3 show derived quantities - the characteristic age, the magnetic field at the light cylin¬ 
der, the spin-down luminosity, and the complexity parameter (oc ~ 5(P/10“^®)^/'^P“®/^"^) as 
presented by Gil & Sendyk (2000). Column 8 gives spectral indices (lUspec) and column 9 
summarizes the X-ray luminosities of the pulsars. These are given for the 2-10 keV band 
unless otherwise stated. 

PSR B1957-I-20 has comparable values of Plc, E, and Oc to the four millisecond pulsars 
that emit giant pulses. Young pulsars like PSR B0540—69 and the Crab have much higher 
values of E and Oc, but they also emit many more giant pulses. If any of these attributes 
dictate giant pulse emissivity, we would expect PSR B1957-I-20 to emit giant pulses. If the 
pulses we see are not giant pulses, then it is plausible that there is another population of 
pulses that has an even lower rate of emission. Presumably these pulses would take the 
form of very narrow spikes that are restricted in pulse phase, just like the pulses we see. As 
invoking two populations of identical looking pulses is contrived, we believe that we have 
seen purely giant pulse emission, or giant pulse emission superimposed on a base of ordinary 
emission. An alternate idea that we do not favor is that at moderate energies ordinary and 
giant pulses are indistinguishable because the two seemingly disparate populations share a 
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common emission mechanism. The ordinary pulses of PSR B1937+21 show no sign of modu¬ 
lation (Jenet & Gil 2004) and the giant pulses only marginally coincide with the envelope of 
ordinary emission. The emission mechanisms are therefore quite distinct for PSR B1937-1-21, 
and consequently it seems unlikely that the population of pulses from PSR B1957-I-20 rep¬ 
resents the transition of a single pulse-population from ordinary-like to giant-like emission. 
At this stage we cannot dehnitely state that PSR B1957-1-20 emits giant pulses. Further 
supporting evidence could be made by establishing power-law statistics and hnding a cor¬ 
relation in phase with an X-ray pulse. The general task of identification of weak pulses as 
giant pulses is more difficult. Evidence could include giant pulses having different DMs to 
ordinary pulses, or characteristic timescales much shorter than those seen for microstructure. 

Although PSR B1957-I-20 has a similar Rlc to the millisecond pulsars that emit giant 
pulses, its emission rate is signihcantly lower. In particular, its rate would appear to be 
~ 100 times lower than PSR B1823—3021A, despite the fact that PSR B1957-I-20 has a 
higher Rlc- Magnetic inclination angle and other geometric factors must play some role, 
but it is difficult to see how they could account for such an enormous difference in emissivity. 
So although the magnetic held at the light cylinder does seem to be a reasonable determinant 
of whether or not a pulsar emits giant pulses, it alone is not a trustworthy indicator of the 
rate of emissivity. 


4.3. Giant Pulse Emitters 

PSR J0218-I-4232 is the fourth millisecond pulsar found that has been shown conclu¬ 
sively to emit giant pulses. All four such millisecond pulsars have high values of E and the 
complexity parameter. The three observed in X-rays are very luminous in the 2-10 keV band 
and have hard photon indices (see Table 3 and references therein). It is tempting to suggest 
that one or more of these characteristics are better indicators of emissivity rates than Rlc- 
However, PSRs B1957-I-20 and J1843—1113 do not have corresponding values that are so 
much lower that E and the complexity parameter can be discriminated from Rlc as the 
primary determinant of whether or not a millisecond pulsar emits giant pulses. In fact Rlc, 
E, and the complexity parameter have such similar P-P dependences that we do not think 
observations of millisecond pulsars can ever discriminate between them. 

Which other millisecond pulsars could emit giant pulses? Pulsars with high Rlc still 
seem to be good candidates, but this parameter no longer appears to guarantee a rate 
sufficiently large to give a high detection count. Table 3 shows that the millisecond pulsars 
that emit giant pulses all have spectral indices much steeper than the average ctspec = —1-9 
spectrum millisecond pulsars found by Toscano et ah (1998). They also have very low 
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characteristic ages and high X-ray luminosities. Perhaps better sources are young or X-ray 
luminous pulsars in globular clusters? Unfortunately the Galactic globular cluster population 
is old and so most cluster pulsars are likely to be too old to be good candidates for giant 
pulse emission. Consider PSR J0024—7204J, which has a 0.5-6 keV X-ray flux of Lx = 

2 X 10^^ergss“^ (Heinke et al. 2005). This is a similar luminosity to PSR B1957-I-20, so we 
do not expect PSR J0024—7204J to emit giant pulses at a high rate. Since PSR J0024—7204J 
has the highest X-ray luminosity of the identihed millisecond pulsars in 47 Tucanae, we do 
not consider 47 Tucanae to be a good candidate cluster for giant pulse emission. The clusters 
most likely to host populations of the young and X-ray luminous millisecond pulsars prone 
to emitting giant pulses are instead those that appear to contain young pulsars, such as the 
core-collapsed clusters M15 and NGC 6624. 

Perhaps all millisecond pulsars emit giant pulses at even lower rates than PSR B1957-I-20? 
The best candidates for verifying this hypothesis are nearby millisecond pulsars that have 
pulses that are not significantly scatter-broadened. Should bright pulsars like PSR J0437—4715 
emit nanosecond-timescale pulses, then high time-resolution studies could potentially probe 
their pulse populations down to very low energies. Such studies could reveal giant pulses 
occuring at rates smaller by factors of ~ 1000 than seen for PSR B1957-I-20. 


5. Conclusions 

We have searched four millisecond pulsars for individual pulses of emission with mi¬ 
crosecond timescales and have found such emission from two of them. Only four individual 
pulses were detected from PSR B1957-I-20 in 8003 s of observations centered at 825 MHz. As 
these pulses are exceptionally narrow there is little scattering-induced pulse-broadening at 
least some orbital phases. Although it is debatable whether or not these strong pulses are 
true “giant pulses”, we can say that the giant pulse emission rate from PSR B1957-I-20 is 
signihcantly less than the rates for other pulsars with similar values of magnetic held at the 
light cylinder. Although Rlc can be used as a rough guide to whether a pulsar emits giant 
pulses, we suggest it is a poor indicator of the emission rate. 

PSR J0218-I-4232 emits giant pulses at a low rate that is inconsistent with the hndings 
of Joshi et al. (2004). It is most likely that the pulses reported by Joshi et al. are spurious. 
The giant pulses of PSR J0218-I-4232 are conhned to two narrow phase regions separated by 
roughly 50% of phase which align in phase with the peaks of the X-ray prohle and roughly 
coincide with the minima of the integrated pulse-prohle in the radio band. This strong 
correlation between X-ray and radio properties conhrms that the two emission processes 
originate in similarly dehned regions of the pulsar magnetosphere. 
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Most of the 139 giant pulses observed from PSR J0218+4232 at a center frequency 
of 857 MHz had relatively low energies, typically only a few times the mean pulse energy. 
Only three had energies above 10(i7) and none had energies above 20(i7). The pulses exhibit 
power-law statistics, are only found in narrow phase windows that coincide in phase with the 
X-ray pulse-components, and are very narrow just like the giant pulses of PSR B1937-I-21; 
it is apparent then that “giant” pulses should be dehned not through large flux densities, 
but by these three properties. The brightest pulse seen at a center frequency of 1373 MHz 
seems to be around 500 ns in duration when viewed at 125 ns time resolution. At higher time 
resolution hner features become apparent, but it is unclear whether these are signihcant. 

PSR J0218-I-4232 is the fourth millisecond pulsar found to emit giant pulses after PSRs 
B1937-I-21, B1821—24, and J1823—3021A. All four have low characteristic ages and steep 
radio spectra. With the exception of PSR J1823—3021A which has not been observed in 
X-rays, the four pulsars all have high X-ray luminosities and exhibit power-law spectra. The 
presence of X-ray emission with a steep power-law spectrum therefore seems to be the best 
indicator of whether a millisecond pulsar emits giant pulses. Radio observations would be 
expected to show that narrow giants will be present at the phase of the X-ray emission. 
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Table 1. Summary of searches for giant pulse emission. 


PSR 

V 

(MHz) 

(MHz) 

^obs 

(s) 

(xlO^) 

{E) 

(Jy ■ /is) 

P/lim 

m) 

A'det 

Notes 

J0218+4232 

857 

128 

3456 

15 

18 

0.58 

24, 75 

(1) 


857 

128 

5216 

22 

9.2 

0.95 

5, 35 

(2) 


1373 

64 

3349 

14 

3.9 

2.5 

4,8 



825 

64 

293 

1.3 

8.7 

1.7 

0,4 


J1012+5307 

825 

64 

720 

1.4 

60 

0.27 

0 



1373 

64 

6077 

12 

15 

0.72 

0 



1437 

64 

6951 

13 

22 

0.51 

0 


J1843-1113 

825 

64 

3198 

17 

6.2 

4.3 

0 



1373 

64 

791 

4.3 

2.3 

5.8 

0 



1437 

64 

1436 

7.8 

1.9 

6.9 

0 


J1959+2048 

825 

64 

8003 

50 

4.7 

4.0 

4 



1373 

64 

4684 

29 

0.72 

16 

0 

(3) 


1437 

64 

4760 

30 

0.63 

19 

0 

(3) 


Note. — (1) 2004 August observation; (2) 2004 October observation; (3) No 
phase-coherent timing solution was available because of incorrect time-tagging. The 
flux density (S') used to derive the given parameters is given by the S = 0.35(z//1490 
MHz)“^ mJy relation of Fruchter et ah (1990). 







Table 2. 


Giant pulse emission rates from a selection of pulsars. 


PSR 

w (MHz) 

Phase range 

K 

a 

P(E > 20{E)) 

Sgp(E > 20(E)) 

Sgp(E > 0.1(E)) 

References 

J0218+4232 

857 

A (August) 

1.2 X 10“® 

1.5 

1.3 X 10“’’ 

8.0 X 10“® 

1.1 X 10“'* 

This work. 

J0218+4232 

857 

B (August) 

5.6 X 10“® 

1.9 

1.9 X 10“’' 

8.0 X 10“® 

9.4 X 10“'* 

This work. 

J0218+4232 

857 

B (October) 

4.4 X 10“® 

1.7 

2.7 X 10“’’ 

1.3 X 10“® 

5.4 X 10“'* 

This work. 

J0218+4232 

1373 

all 

3.4 X 10“® 

1.5 

3.8 X 10“’’ 

2.3 X 10“® 

3.2 X 10“'* 

This work. 

81957+20®- 

825 

main pulse 

2 X 10“'' 

2 

5 X 10“® 

2 X 10“® 

4 X 10“^ 

This work. 

60531+21 (Crab) 

146 

main pulse 

2.8 X 10“1 

2.5 

1.5 X 10“** 

5.1 X 10“® 

n/a 

(1) 

B0531+21 (Crab) 

800 

all 

9.8 

2.4 

8.3 X 10“3 

2.9 X 10“' 

n/a 

(2) 

B0540-69 

1390 

early 

2.4 X 10“'' 

1.5 

2.7 X 10“** 

1.6 X 10“'' 

n/a 

(3) 

B0540-69 

1390 

late 

7.6 X 10“' 

2.1 

1.4 X 10“^ 

5.3 X 10“'' 

n/a 

(3) 

B0540-69 

1390 

all 

2.6 X 10“' 

1.8 

1.2 X 10“^ 

5.4 X 10“'' 

n/a 

(3) 

B1937+21 

430 

all 

3.2 X 10“'' 

1.8 

1.5 X 10“** 

6.6 X 10“® 

n/a 

(4) 

B1937+21 

1650 

all 

2.8 X 10“'* 

1.4 

4.2 X 10“® 

3.0 X 10“** 

2.5 X 10“® 

(5) 


^Rates are indicative only due to the very small number of pulses analysed. 


References. — (1) Argyle Gower (1972); (2) Lundgren et al. (1995); (3) Johnston et al. (2004); (4) Cognard et al. (1996); (5) Soglasnov et al. 
(2004). 







Table 3. Pulsar characteristics. 


PSR 

P 

P 

T 

^LC 

E 

<3.0 

•^spec 

^X^-10 keV) 

References 


(ms) 

(10"^“-) 

(Myr) 

(lO^G) 

(10^^ ergs s~^) 



(10^^ ergs s~^) 


J1823-3021A® 

5.44 

3390 

26 

25 

810 

28 

-2.7 

Unknown 

(1,2) 

B1821-24® 

3.05 

1620 

30 

74 

2200 

33 

-2.3 

13, 12.8 (0.5-8keV) 

(1,3,4,5) 

B1937+21 

1.56 

106 

230 

102 

1100 

23 

-2.6 

0.5-5.7 

(1,3,6,7,8) 

J0218+4232 

2.32 

77.4 

480 

32 

240 

16 

-3.0 

1.3 (1-lOkeV); 1.2-1.6 

(1,9,10,11,12) 

J1012+5307 

5.26 

9.73 

8600 

1.5 

2.6 

6.3 

-1.9 

0.003 (0.2-10keV) 

(13,14,15) 

J1843-1113 

1.85 

9.59 

3100 

20 

59 

11 

Unknown 

Unknown 

(1) 

B1957+20 

1.61 

11.5 

2200 

31 

110 

14 

-3.0 

0.16 (0.5-7keV) 

(1,6,16,17) 


^'Parameters ignore acceleration in the gravitational potential of the host cluster. 


References. — (1) Manchester et al. (2005); (2) Toscano et al. (1998); (3) Foster et al. (1991); (4) Becker et al. (2003); (5) Mineo et al. 
(2004); (6) Toscano et al. (1999); (7) Takahashi et al. (2001); (8) Nicastro et al. (2004); (9) Navarro et al. (1995); (10) Mineo et al. (2000); 
(11) Kuiper et al. (2002); (12) Webb et al. (2004a); (13) Nicastro et al. (1995); (14) Lange et al. (2001); (15) Webb et al. (2004b); (16) 
Fruchter et al. (1990); (17) Stappers et al. (2003). 
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Fig. 1.— Top: Phases and energies of pulses detected from PSR J0218+4232 in the 2004 
August observation are superimposed on an integrated pulse prohle. Bottom: The Chandra 
HRC-S 0.08-10 keV pulse prohle of PSR J0218-I-4232 (Kuiper et ah 2004) has been phase- 
aligned with the radio prohle using the absolute timing of Rutledge et al. (2004). 
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Fig. 2.— Cumulative distribution of giant pulse energies for observations of PSR J0218+4232 
centered at 857 MHz when viewed in terms of absolute energy (left panel) and relative to 
the mean pulse energy (right panel). 
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Fig. 3.— The intensity of the central portion of the strongest giant pulses from 
PSR J0218+4232 and PSR B1957+20 when seen with a sampling interval of 15.625 ns (solid 
line) and 125ns (raised dotted line). 
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Pulse phase 


Fig. 4.— Relative phases of the individual pulses detected from PSR B1957+20. The top 
four panels show the intensities of the pulses and the bottom panel shows the intensity of 
the integrated pulse prohle. 




